Ultra-sensitive Diamond Magnetometry Using Optimal Dynamic Decoupling 
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New magnetometry techniques based on Nitrogen Vacancy (NV) defects in diamond have received 
much attention of late as a means to probe nanoscale magnetic environments. The sensitivity of a 
single NV magnetometer is primarily determined by the transverse spin relaxation time, T2. Current 
approaches to improving the sensitivity employ crystals with a high NV density at the cost of spatial 
resolution, or extend T2 via the manufacture of novel isotopically pure diamond crystals. We adopt a 
complementary approach, in which optimal dynamic decoupling techniques extend coherence times 
out to the self-correlation time of the spin bath. This suggests single spin, room temperature 
magnetometer sensitivities as low as 5pTHz~^''^ with current technology. 
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Room temperature single spin magnetometry using the 
nitrogen-vacancy (NV) centre in diamond has the poten- 
tial to revolutionize nanoscale imaging through funda- 
mentally new detection modes [H-Q. Proposals to image 
nanoscale environments exhibiting static (DC) and oscil- 
latory (AC) magnetic fields using NV systems 0, 01 have 
since been demonstrated experimentally By ex- 

ploiting measured changes in quantum decoherence 0, II] 
these techniques have been extended to include more gen- 
eral classes of randomly fluctuating (FC) magnetic fields, 
with comparable sensitivity to the AC case j8|- One 
prominent reason such room temperature single spin de- 
tection techniques are of interest is because of the po- 
tential to develop fundamentally new imaging modes for 
biological systems with nanometer resolution Q. 

The sensitivity of an NV magnetometer is governed 
by the transverse spin relaxation (dephasing) time T2, 
which for AC detection using isotopically-pure diamond 
has been demonstrated at 4nTHz~^^^ jfl]. Many im- 
portant detection problems, particularly those related to 
biology stand to gain significant improvements with 
increased coherence times. In this paper we show how 
optimal dynamic decoupling techniques can be ex- 
ploited to increase sensitivities by over two orders of mag- 
nitude. With sensitivities in the pTHz""'^''^ regime, and 
nanoscale spatial resolution, the ultra-sensitive NV mag- 
netometer proposed here may have profound implications 
for nano-bio imaging and sensing. 

It is well known that coherence times may be improved 
with the use of concatenated (CDD) lllll. random and pe- 
riodic dynamic decoupling schemes [ll|, UM,. The more 
recent Uhrig (UDD) scheme was shown to be optimal 
for decoupling a spin qubit from a bosonic bath [10| , and 
has since been applied to a much broader class of sys- 
tems in which dephasing is the dominant decoherence 
chanriel, including classical noise and spin bath sys- 
tems [IJ] . This optimality stems from the fact that the 
required resources of the UDD sequence grows linearly 
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FIG. 1: (color online) (a) NV-centre diamond lattice defect, 
(b) NV spin detection through optical excitation and emission 
cycle. Magnetic sublevels nis = and rus — ±1 are split by a 
D=2.88 GHz crystal field. Degeneracy between the rUs — ±1 
sublevels is lifted by a Zeeman shift, Suj. Application of 532 
nm green light induces a spin-dependent photoluminescence 
and pumping into the — ground state, (c) Examples of 
controlled- AC fields (solid) as seen by the NV centre (dashed) 
in the presence of the 1"* , 2°"^ and 5**^ UDD sequence. Nega- 
tive regions of the AC trace are mapped to positive, ensuring 
maximal phase accumulation of the NV spin. Slow FC fields, 
such as the surrounding nuclear spin bath, will be suppressed, 
permitting the detection of AC and fast FC fields with greater 
sensitivity. 



with the order to which the environmental effects are 
suppressed [141 • However, UDD is unable to suppress 



longitudinal relaxation, whereas CDD can, but does so at 
exponential cost. The use of UDD is ideally suited to the 
NV centre owing to long relaxation times (Ti > Is [Gi]), 
even at room temperature. The large Debye tempera- 
ture allows for negligible decoherence via phonon excita- 
tion of the crystal lattice, and a large zero field splitting 
(2.88 GHz) of the ground state magnetic sublevels pre- 
vents longitudinal spin-spin relaxation. Hence, for the 
NV centre, longitudinal relaxation may be neglected, en- 
suring UDD is the optimal decoupling method. 

Both AC and FC magnetometry schemes are based 
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upon a spin-echo microwave control sequence, in which a 
TT pulse is used to flip the qubit at the half way point of 
its evolution (Fig.[T][c)), suppressing any quasi-static ef- 
fects of the spin bath. The AC scheme is concerned with 
detection of fields of the form 5ac = bo sm{i't), where the 
TT pulse coincides with t ~ tt/i/, ensuring a non-zero inte- 
gral of the field trace, and hence maximal phase shift of 
the NV spin. In this paper, we incorporate UDD into the 
NV magnetometry scheme and show that the sensitivity 
of a single NV centre magnetometer may be as low as 5 
pTHz^^/^ with the use of existing technology. 

Recently it was shown that a particularly intuitive 
analysis of a spin qubit placed in a slowly fluctuating 
classical magnetic field could be performed by expanding 
the time dependent field as a Taylor series [8j. Whilst 
this may seem like a special case, this technique applies 
to a more general class of pure-dephasing quantum prob- 
lems in which longitudinal relaxation may be ignored. 
We investigate the effect of UDD on a spin qubit placed 
in such a field, and show that it is the n"^ UDD sequence 
that suppresses the effect of all terms up to and includ- 
ing order n in the Taylor expansion of the field. These 
results are applied to the NV centre and used to obtain 
improved sensitivities for NV based magnetometry. 

An NV centre aligned along the z axis interacting with 
an external magnetic field is described by the Hamilto- 
nian H — Hzis + Hext + "Hint- The first term describes 
the zero field splitting of the ground state Zeeman levels, 
-H^f, = hDS'^, where £) = 2.88 GHz. The interaction with 
an external magnetic field Boxt(i) is described by 'Hcxt- 
The fields we consider here are small relative to D and are 
hence unable to induce a spin-flip, permitting us to ignore 

all Sx,y terms, giving Hext = 7nvS • Bext « Tnv^oxt^'z, 
where 7 is the electron gyromagnetic ratio. For simplic- 

(z) 

ity, we put Sg^t = ^oxt- The flnal term describes the 
interaction with the paramagnetic environment of the 
diamond crystal. As we will see, for resolving the non- 
unitary dynamics of the reduced density matrix of the 
NV centre, these interactions may be subsumed into a 
single 'internal magnetic field,' B[^t{t). We define the 
fluctuation regime of the external/internal environment 
via the dimensionless numbers, Ooxt = (7fo'^*''cxt)~"'^ and 
Bint = (7crj,"*Ti„t)~\ where crg'^Vcrj,"* are the RMS held 
strengths, and Text /Tint are the correlation times of the 
internal and external environments. Rapidly and slowly 
fluctuating flelds satisfy 6^1 and 6^1 respectively. 

An arbitrary time-dependent magnetic field may be 
decomposed as a Taylor series in t: B{t) = Y^^=o^kt^- 
The validity of this expansion rests upon the condition 
that a„_)_ii"+^ < a„i" for all n. As shown in 8] this is 
satisfied for t < T-mt/V^- For times t ^ Tint the qubit 
will exhibit motional-narrowing behaviour. 

In many cases of practical interest, Bint is the sum 
of fields from a large number of dipoles. This im- 
plies that the {aj} are normally distributed, with zero 



mean at room temperature and variance cr^ 



This leads to the following dephasing envelope, V(t) 
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Since Fq ^ F/cVfe > 1, Fq serves to define the free 
induction decay time, = I/Fq. For a 1.1% "C 
bath, the variance of the magnetic field is given by 
ctq — ^.i{Bf). Approximating this sum by an integral 

gives (To « y^l^nc^cMN ~ 2^T, and Tj* « 4/^s, in 
good agreement with 0, [l^. Similarly, using an 0.3% 
^^C bath yields T| « 15 /its, in agreement with 

The Hahn-echo sequence removes the effect of a static 
field on the system (though there are an infinite number 
of sequences that achieve this) , as each tt pulse effectively 
sends B — )• —B. For a field described by J2k^kt'^' 
Hahn-echo sequence will remove the effect of the on term. 



and modify all other terms as 1— (1 — 2 



For 



an NV centre, the correlation time of the environment is 
dictated by interactions between ^"^C nuclei. A straight- 
forward calculation shows Tint ^ \f^^f^ /(^(^9cI^n) — 
15 ms. Substituting this into Eq. ([T]) for j — 1 gives 
Fi = 2.1 kHz. We identify T2 = 1/Fi = 400 /xs, again 
in agreement with experiment 0, [S]- For an 0.3% ^^C 
bath we achieve T2 = 1.5 ms as seen in Q. We do not 
define Tint via interaction between ^'^C nuclei and any 
background fields, i?o, since this manifests as decays and 
revivals on timescales of t^ ^ 7c and does not repre- 
sent a true loss of information. 

We may wish to modify our § — pulse sequence 

in order to remove the effect of the a\t term. If we ap- 
ply pulses at t = J and t = we find that the ef- 
fects of both qq and Oi terms are suppressed. In gen- 
eral, suppression of all terms up to and including or- 
der n will require n + 1 7r-pulses. We define Tn^k as the 
time at which the k^^ pulse is applied in the sequence 
that suppresses all field components up to, and includ- 
ing, order n. Determination of the n -I- 1 elements of 
the set Vn — {t„,o , • ■ • , Tn,n+i} will require the solu- 
tion of the following set of n -I- 1 algebraic equations for 
T„.fc. For ao, 2 (-1)'"' r„,k + (-1)"+' r = 0; for 

«i, 2 Eltl i-l)'-' rl, + i-ir' = 0; up to 



fe=l 



n+1 _ 



= 0, (2) 



To this point, there is some freedom in our choice of pulse 
sequence, as Eqs. ([2]) are satisfied to n = 1 by both CDD 
and UDD. However, if we solve Eqs. ([2]) to n = 2, wc find 
V2 = If (2 - 72) , §, § (2 + V2)}, which is the third 
UDD sequence. It is reasonable to conjecture that the 
n^^ component of the field, a„i", will be suppressed by 
the {n + UDD sequence, which we prove below. 
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For an interrogation time of r, the time of application 
of the fcth pulse in the nth UDD sequence is given by 
T„,fc = Tsin^ (2^) M, where I < k < n. We wish 

to show that the (n + 1)"^ UDD sequence suppresses the 
effect of terms up to and including the nth term in the 
Taylor expansion of a time dependent magnetic field; and 
that the effect of all terms beyond n will be reduced. 

The phase accumulation of a spin qubit is propor- 
tional to the time integral of the magnetic field to which 
it is exposed. Each tt pulse exchanges the basis states 
of the qubit Hilbert space, which has the same effect 
on their relative phase as mapping B M- — _B. Re- 
call from Eq. ^ that the action of the (n + 1)*^ Uhrig 



pulse sequence 



will be 



to modify each of the a, via 
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To prove our claim, we must firstly show that. 
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(-1)" 



(3) 



Using N = n-|-2, sin(x) = ^ (e*^ — e and expanding 
as a binomial series in j, we arrive at 
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where ajmN — "^aT^ ■ Note that we have added the 
fc = term, since X]j=o (^1") i^^V — ^- Since we have 
restricted ourselves to to < A^— 1, we have that e'^^^J™" 7^ 
1, so we are free to sum over A: as a geometric series in 
Eq. (U), which is not possible for m > N . This gives 

of terms in the (2to -I- 1*'') row of Pascal's triangle, and 
evaluates to 2^™, proving Eq. ([3|). 

Replacing N with n + 2, we can see that i-> , V j < 
n. Hence all terms in the Taylor expansion of B(t) up 
to and including order n are zero. Furthermore, since 
sin^(a:) < 1, the term in the square brackets is always less 
than 1, hence the effect of the remaining aj is reduced. 

Two immediate consequences are that the (n -I- 1)*'' 
UDD sequence will suppress dephasing to n"^ order; and 
that dephasing effects beyond n*^ order will be reduced. 
If ao, . . . , a„ = 0, then (ag) , . . . , (al) = 0. By Eq. © 
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FIG. 2: (colour online) (a) Dephasing envelopes, V^"\ for an 
NV centre in a 1.1% ^^C bath for n tt pulses. As n — >■ 00, 
approaches the Heaviside step function, H{rcxt —t). (b) 
Effect of the number of tt pulses on NV coherence times for 
different ^^C concentrations. In each case, the coherence time 
is limited by Tint- 



with the UDD sequence, or an AC field of controllable 
frequency whose nodes coincide with each tt pulse, which 
could be realised by a single spin or ensemble of spins be- 
ing driven by a controllable microwave field. Sensitivity 
to rapidly fiuctuating fields will also be improved, since 
fields with correlation times shorter than the interroga- 
tion time will not be refocussed by the UDD sequence. 

We denote the dephasing envelope in the pres- 
ence of the (n + l)th pulse sequence as V'^^^t) = 

nfc^n^i"^(0 [Fig.[2lja)]. In the presence of background 
dephasing described by Eq. ([IJ , the minimum induced 
phase from i3ext(i) that may be measured is A0(&) = 
[C\/N'D^"\t)]~^ , where C describes photon shot noise 
and imperfect collection Q, and N is the number of 
measurements taken. Typically C < 0.3, however vast 
improvements have recently been demonstrated by en- 
tangling the NV spin with proximate nuclear spins, per- 
mitting repetitive readout of the NV spin state [3: [13 ■ 
We now discuss the relevant detection protocols and 
associated sensitivities for different fields to which these 
techniques apply. For a telegraph signal (square wave- 
form) switching between ±i?o in sync with each tt pulse, 
the qubit will acquire the maximum possible phase 
for a given interrogation time, A0 = 'jBqt. This 



T 



gives a magnetic field sensitivity of rjl^ = Bq 
[C7v^X>(")(t)]~^ For all cases where Oint < 



r„, are zero. In what follows, we ex- ^^""^ t^^^* ^fe > Ffc+i, and Ffc > F 



all rates, Fq, . 

amine how these ideas may be implemented to improve 
the sensitivity of a qubit magnetometer. 

From the above analysis, we see that the application 
of UDD sequence of any order will decrease the NV de- 
phasing rate due to the internal fluctuating magnetic 
field. This allows us to extend the interrogation time, 
and hence improve the sensitivity to an external magnetic 
field, Bext(t). Clearly the dynamics of Bcxt will be an im- 
portant factor in ensuring that the effect of the external 
field is not also suppressed by the pulse sequence. Sim- 
ple examples include telegraph signals switching in sync 
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may approximate the total dephasing envelope, P*^"^ , by 
the its leading order contribution. That is 2?("^(i) ^ 



exp 



(F„+ii) 



2n+4 



time is T = F^^jjij^ (4n 



, implying the optimal interrogation 
ji/(2n+4)^ We then find the sen- 
sitivity to be bounded above by j]'^''^ < ^\//e0inP"- 
This upper bound, together with the actual sensitivity 
(see below), is plotted in Fig.tSja) for an NV. Notice that, 
as 9int — 1, there is little to be gained by applying UDD. 

We now compute the sensitivity of the probe by taking 
into account the effect of the n**^ order pulse sequence on 
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FIG. 3: (colour online) (a) Fundamental shot-noise sensitivity 
limits assuming C = 1 and a synchronised square wave-form 
for different ^'^C concentrations. Sensitivities to controlled AC 
signals are modified by a factor of 7r/2; those due to FC signals 
satisfying 0cxt > 1 are modified by a factor of 20oxt • Smooth 
curves show the analytic upper bound on the sensitivity, r;^*' . 
(b) Effect of finite pulse widths of 50 ns on sensitivities. The 
optimal sensitivity occurs for n = 13 tt pulses. 



the n}^ order Taylor coefficient. This results in improved 
sensitivity beyond that indicated above, as a reduction 
in the Uk leads to a reduction in the total decoherence 
rate, and hence an extended interrogation time. 

The dephasing rates are found via Eq. ([1]), where 



An) 



2ES(-1)* 



^isin2^+2 



J, will be given by the solution to 

2k+2 

— 1 = for i, and are plotted against 



The actual dephasing time due to the combined ef- 
fect of all the r^"' 

l^k=n+l k ''J 

the number of pulses used in the sequence in Fig.[21[a). 
From this we see that dephasing times asymptote to the 
correlation time of the bath. This is to be expected, since 
information lost to the bath cannot be recovered by con- 
trol of the NV spin alone. Retaining all terms, the sensi- 
tivity is then yy^"^ (r) = ^i^^ exp 



2k+2 



By minimising this function with respect to r we deter- 
mine the optimal sensitivity, ryt"^ shown in Fig.|31Ja). 

In many proposals magnetic resonance techniques 
are used to drive the sample magnetisation at some con- 
trolled frequency. By synchronising frequencies with the 
chosen pulse sequence, a piecewise continuous sinusoidal 
signal may be produced [Fig.[Tfc)]. The sensitivity is 
hence modified by a factor of 7r/2, ??ac(''') = f ??t"^- 

Instead of sensing a controlled signal, we may wish to 
characterise a random, rapidly fluctuating (0ext ^ 1) ex- 
ternal field, -Bext ■ Rather than yielding an identical phase 
shift in each run, such a field will introduce an additional 
dephasing component and may be detected as a pertur- 
bation in the dephasing rate [8| . The dephasing envelope 
will be modified by a factor of I?cxt {t) = exp (— Foxti), 
where Foxt = 57pCext''ext- The sensitivity with which 

(Text = \J (-^oxt) ~ (5cxt)^ niay be measured is then 

VSii''') — 26cxt'7tr\8], making the field more difficult to 
detect as the fluctuation rate increases. This is consistent 



with motional narrowing phenomena in NMR, in which 
high frequency noise is known to have a reduced effect 
on the sample T2 as compared with quasi-static noise. 

Coherent manipulation of the NV is achieved via 
a Rabi cycle, thus instantaneous tt pulses cannot be 
achieved in practice, and lead to additional decoher- 
ence effects. For a Rabi frequency of il, the decoher- 
ence envelope is given by = [1 -I- (7^0-0 
[l8| , and typical pulse errors are ~ 1% Hence for 
n TT pulses, the sensitivity will be worsened by a fac- 
tor of A?7 ~ 0.99"+i[l + ^{y/^-fao/n)% as shown in 
Fig-EIb). For a pulse width of 50 ns, 13 tt pulses is found 
to be optimal, with T]il^'^ « 5.5pTHz-i/2 

We have theoretically investigated the improvements 
associated with the application of the optimal Uhrig con- 
trol sequence to an NV-based magnetometer. Results 
show that dephasing times due to the paramagnetic im- 
purities are ultimately limited by the self correlation time 
of the fluctuating environment, thus NV magnetometer 
interrogation times may be extended by nearly four or- 
ders of magnitude beyond the free-induction decay time. 
In light of these results, we have shown that incorpora- 
tion of UDD into current single NV magnetometer pro- 
tocols may yield sensitivities below 5pTHz~^/^ at room 
temperature in the near future. Such techniques have 
the potential to yield great improvements to nano-scale 
sensing, particularly nano-biological processes occurring 
at room temperature. 
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